On the dimerized phase in the cross-coupled antiferromagnetic spin ladder 
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We revisit the phase diagram of the frustrated s=l/2 spin ladder with antiferromagnetic rung 
and diagonal couplings. In particular, we reexamine the evidence for the columnar dimer phase, 
which has been predicted from analytic treatment of the model and has been claimed to be found 
in numerical calculations. By considering longer chains and by keeping more states than in previous 
work using the density-matrix renormalization group, we show that the numerical evidence presented 
previously for the existence of the dimerized phase is not unambiguous in view of the present more 
careful analysis. While we cannot completely rule out the possibility of a dimerized phase in the 
cross-coupled ladder, we do set limits on the maximum possible value of the dimer order parameter 
that are much smaller than those found previously. 

PACS numbers: PACS number: 75.10.Jm 



I. INTRODUCTION 

In spin chains, frustration due to competing an- 
tiferromagnetic nearest- and comparably strong next- 
nearest-neighbor coupling leads to the appearance of a 
new gapped, dimerized phase^'^ whose properties are 
markedly different from that of the uniform spin-liquid 
state of the antiferromagnetic spin-1/2 chain. A much 
richer phase diagram is found in spin ladders, where the 
couplings between nearest-neighbor spins situated along 
the legs and on the rungs of the ladder may compete. 
These phases and the quantum phase transitions between 
them have been extensively studied recently, both exper- 
imentally and theoretically."^ The theoretical approaches 
include both analytic and numerical work. The analytic 
treatment of the bosonized form of the spin Hamiltonian 
has been widely used to predict the phase diagrami^"— 
There is general agreement that, in two-leg spin ladders, 
the most relevant operator generates either a topologi- 
cally ordered rung-singlet or a Haldane phase, depend- 
ing on the relative strength of the intraleg and interleg 
couplings. 

The competition between the couplings is, however, a 
delicate problem when the intraleg and the interleg cou- 
plings are of the same order of magnitude, and the type 
of order is, in some cases, determined by marginally rel- 
evant terms generated within the renormalization group 
(RG) procedure. Among the many terms generated, op- 
erators leading to dimerization appear. On this basisj^ 
it was proposed that a dimerized phase appears between 
the topologically ordered rung-singlet and the Haldane 
phases for the so-called cross-coupled ladder, in which 
spins are coupled not only along the legs and across the 
rungs but also diagonally, i.e., on opposite legs and on 
different, neighboring rungs. 

The existence of an analogous dimerized phase induced 
by four-spin exchange terms has been proposed some 
time ago.— This dimerized phase has indeed been found 
to exist numerically^^"— It is also well-established that 
a dimerized phase exists when a second-neighbor inter- 



action on the legs of the ladder are present^! or when 
the interleg coupling is ferromagnetic)^ The existence of 
the dimerized phase has, however, been disputed for lad- 
ders with pairwise spin interactions orAj^— while other 
work has presented evidence supporting the presence of 
this phase. ?2i2g, 

The density-matrix renormalization group (DMRG) 
procedure invented by White^'' is a particularly use- 
ful numerical method to study such problems because 
it makes possible the high-precision calculation of the 
ground-state properties and low-lying excitations of low- 
dimensional quantum systems. The DMRG has indeed 
been widely use d^^'^^" — to explore the possible phases 
of spin ladders. The length of the ladders from which 
finite-size scaling is done and the accuracy of the DMRG 
calculation can be, however, decisive. Motivated mainly 
by the recent work of Hikihara and Starykh^^, we have 
undertaken a careful analysis of the numerical data to 
see if there is indeed compelling evidence for the dimer- 
ized phase in the cross-coupled ladder using the DMRG 
procedure. 

The paper is organized as follows. In Sec. II, we 
summarize earlier results on the phase diagram of the 
cross-coupled frustrated antiferromagnetic spin-1/2 lad- 
der. The aspects of the DMRG procedure that deter- 
mine its accuracy are briefiy discussed in Sec. III. Our 
numerical results are presented in Sec. IV. Finally Sec. V 
contains our conclusions. 



II. PREVIOUS RESULTS ON THE 
CROSS-COUPLED SPIN LADDER 

We treat a two- leg spin-1/2 ladder with pairwise 
Heisenberg couplings along the legs, across the rungs, 
and diagonally between sites on adjacent rungs, but on 
opposite legs. The Hamiltonian can be written 
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where 

T-liS+l — J\\ {Sl/i ■ Sl^i+l + S2,i • 82,1+1) 

+ ^J± {Si,i ■ 82.1 + 8i,i+i ■ 82,1+1) (2) 

+ Jy. (8i,i ■ 82,1+1 + 82,1 ■ 8i,i+i) , 

and 8i,i and 82.1 are the spin operators on rung i and legs 
1 and 2, respectively. A depiction of the couplings of the 
cross-coupled ladder is shown in Fig. [1] Here we will con- 




FIG. 1: (Color online) Depiction of the cross-coupled ladder. 

sider only antiferromagnetic coupling and set the intraleg 
coupling to unity, J|| = 1 which we take as the energy 
scale in the problem. For Jx = 0, it has been established 
that the entire half-line < Jj_ < 00 of the ground- 
state phase diagram, is continuously connectedj^i^ i.e., 
no phase transition occurs for any finite value of J+. 
Since the configuration that dominates the ground-state 
wave function is the product of rung singlets, this phase 
is often referred to as the rung-singlet (RS) phase. Con- 
versely, for J_L = 0, the entire half-line < Jx < 00 
is continuously connected to the Haldane phaseS^ of the 
spin-1 chain and the system is known to have a gapped 
spectrum of spin-1 magnonsj^^- 

When both J± and Jx are finite, frustration due to 
competition between the intraleg and the rung or diag- 
onal interleg couplings potentially leads to the forma- 
tion of new phases. Indications for such new phases are 
given by analytic calculations which can be performed 
for weak interchain couplings using the bosonized form 
of the Hamiltonian. The leading term in the Hamiltonian 
is the coupling between the staggered magnetizations of 
the legs: 

(J_L - 2Jx)ni • n2 ■ (3) 

This term generates a gap in the spectrum irrespective 
of the sign of the coefficient J± — 2Jy . Rung singlets are 
generated for J^ — 2Jx > 0, while the Haldane phase 
appears for J+ — 2 Jx < 0. The first numerical results on 
the model^ii^ were consistent with a first-order transition 
for both weak and strong interchain couplings. Subse- 
quent work has suggested that the transition is actually 
continuous when the interchain coupling is weak, becom- 
ing first-order only for stronger interchain couplingi ^^i^^ 
This scenario has been supported further by an analysis 
based on quantum information entropies. Although two 
critical values of the coupling were found in the two-site 
entropy function for finite systems in this work, it was 
shown that they scale to a single value in the thermody- 
namic limit. It was therefore argued that a direct tran- 
sition between the rung singlet and the Haldane phases 
takes place. 



When J± = 2 Jx , the most relevant term is suppressed 
at the transition between the two gapped phases and 
other, less relevant terms may play a role. In particu- 
lar, such terms lead to the prediction of a spontaneously 
dimerized phasc j^^'^^i^^ The boundaries of this dimerized 
phase have been estimated to lie at 

J± = 2Jx - ^ Jx and J+ = 2Jx - ^ J^ . (4) 

The most recent numerical results, however, indicate a 
much narrower range for the existence of the dimerized 
phase In Ref. [2y the dimerized phase was claimed to 
exist in the range 0.373 < J+ < 0.386 for Jx 0.2. The 
authors of Ref. 22 chose a single point in the middle of 
this range, namely J+ = 0.38, Jx ~ 0.2 to demonstrate 
the existence of the dimerized phase. In what follows, 
we will also restrict ourselves to the vicinity of this single 
point in parameter space and will carry out a careful 
numerical analysis of the dimer order parameter in the 
range 0.36 < Jj_ < 0.4 for Jx = 0.2. 

III. NUMERICAL PROCEDURE 

A. Factors determining the accuracy of the DMRG 
procedure 

In the standard DMRG procedure, the chain is divided 
into two blocks and two extra sites. As the two sites are 
incorporated into larger blocks during the steps of the 
rcnormalization procedure, only a fraction of the new 
block states are kept; the others are discarded. Usually 
the number of block states, M, kept in the truncation 
procedure ranges from a few hundred to a few thousand 
and is set according to the accuracy required and the 
available computational resources.—"— 

As we have pointed out in previous work, keeping the 
number of block states constant for increasingly large 
system sizes leads to increasing error f^ML i.e., the scal- 
ing of block entropy with system size is not taken into 
account In a consistent calculation, the number of 
block states must be increased as system size is increased. 
The finite-size scaling to infinite system size must be 
accompained by a systematic increase of M, and the 
M 00 limit must be taken. We have shown that the 
proper choice of M can be determined by the so-called 
Dynamic Block State Selection (DBSS) approachjiMi in 
which the threshold value of the quantum information 
loss X is fixed a priori. In this approach, the number 
of block states to be kept is obtained from the criterion 
s — sTr""'= < -)^^ where s and sTr""'= stand for the block en- 
tropy before and after the truncation, respectively. The 
minimum number of block states M^-m is set prior to the 
calculation, and the maximum number of block states 
needed to achieve the desired accuracy, Mmax, is moni- 
tored during the DMRG process. This allows for a rigor- 
ous control of the numerical accuracy and gives a stable 
extrapolation to x = 0. 
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B. Choice of a proper order parameter: advantages 
and limitations 



In an earlier worki^ we used the generalization of the 
so-called hidden topological order^ O, to two-leg lad- 
ders, 
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hm ({Sl, + Sl^) 



(5) 



X exp 
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{Sji + Sj2) 



and 



Oevcn = - hm (5f+i 1 + 5f 2) 
\i-j\->-co \ 



(6) 



X exp 



ITT (^f+1,1 + 3^2) 

l=i+l 



to identify various spin-liquid phases.'^* The advantage 
of these specially constructed order parameters is that 
only one of them can be finite in a given phase and they 
both vanish for critical systems. Unfortunately, the rung- 
singlet and dimerized phases would fall into the same 
topological sector because Ocvcn is finite for both cases. 
Thus, they cannot be distinguished by the topological 
order. 

A similar problem arises for the so-called z 
operators which can also be used to distinguish 
different odd- or even-parity valence-bond-solid (VBS) 
states. The operator 



^rung 



exp 



(7) 



defined for ladders with periodic boundary conditions, 
converges to in the rung-singlet phase and to —1 in 
the Haldane phase, while the operator 



exp 



.27r 



N 



E^(^iVi+^2,/) 



(8) 



converges to —1 and +1, respectively, in the two topo- 
logically distinct phases. They give a clear indication 
of the parity of the number of singlet bonds "cut" by a 
line between adjacent rungs, but they do not distinguish 
between the columnar-dimer phase and the rung-singlet 
phase. 

Since we are looking for a possible columnar-dimer 
phase only, a more natural quantity to examine is the 
local dimer order parameter {D) defined as 



D= lim |i:'(iV)l 



(9) 



with 



D{N) = {nN/2-l,N/2) - {'HN/2.N/2+1) ■ (10) 



Thus, dimerization is measured as the alternation in the 
bond energy in the middle of a sufficiently long open 
chain. Here D is finite when the translational symmetry 
of the Hamiltonian is spontaneously broken in the ther- 
modynamic limit. It is worth noting that our definition 
is slightly different from the one used in Ref. [13 for the 
columnar phase, where only the intraleg part of the spin 
couplings was included. Although, quantitatively the 
two definitions give different values for the order param- 
eter, the qualitative behavior is the same. 

For noncritical models with a finite correlation length, 
the end effects decay exponentially and the local quantity 
D{N) is expected to vary to leading order as 



D{N) = D + aN~^ exp(-7V/2^). 



(11) 



This behavior is qualitatively similar to the scaling of the 
gap for periodic boundary conditions (PBC), except that 
the scaling variable is the distance of the middle of the 
chain from the boundary, N/2, and the exponent /? of the 
algebraic prefactor is, a priori, unknown. Here ^ stands 
for the correlation length, which is finite for gapped sys- 
tems. The exponential convergence of local quantities 
such as the dimer order parameter D{N) makes the ex- 
trapolation to the thermodynamic limit more reliable4^ 
In cases when the available data for D{N) not show a 
tendency to go to a finite value as a function of l/N, 
we also carry out extrapolations without the exponential 
term in Eq. (jlip . In addition, in order to obtain an upper 
bound for D, we also use a second order polynomial fit 
in Since the numerical accuracy is of crucial im- 

portance in the present study, we analyze the behavior 
of the dimer order parameter as a function of M and x 
systematically. 



IV. NUMERICAL RESULTS 

As mentioned above, the dimerized phase exists, if it 
exists at all, in a narrow region around J± — 0.38, Jx — 
0.2. We first calculate the length-dependent D{N) using 
open boundary conditions (OBC) on systems of up to 
TV = 512 rungs in steps of 32. We perform four DMRG 
sweeps with fixed numbers of block states, M = 256 and 
512. Note that the longest system studied in Ref. [l^ 
had N — 192 rungs. The results for the dimer order 
parameter at J± = 0.38, Jx = 0.2 are shown on a log-log 
scale in Fi g, p to allow for direct comparison with Fig. 
10 of Ref. 122- The data points calculated at fixed M, 
M — 256 or M — 512, exhibit a clear tendency to go to a 
finite value as a function of N in the infinite-chains limit. 
The flat region develops for M = 256 at chain length 
(N ~ 128) and for M = 512 when N exceeds about 
256. This flattening of the dimer order parameter with 
increasing N was interpreted as evidence for the presence 
of a dimerized phase. 

In order to see what happens for even larger values of 
M, we have carried out calculations keeping M = 1024 
block states for ladders of up to = 256 rungs. As 
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FIG. 2: (Color online) Dimer order parameter D{N) plotted 
as a function of on a log-log scale for various fixed numbers 
of block states (M) for J± = 0.38, Jx = 0.2. The dotted 
line indicates — 192, the largest system studied previously. 
The M — oo extrapolation is taken from the fits depicted in 
Fig. Ha). 

can be seen in Fig. [2l the dimer order parameter does 
not tend to a finite value with N at the available system 
sizes. To carry out the standard extrapolation scheme, 
we have fitted the M-dependent D{N) for various system 
sizes as a function of 1/M using the function 

D{M, N) = D{N) + a{N)A'r^'^^'^ , (12) 

where D{N), a{N), and f3{N) are free fitting parameters. 
These fits are shown in Fig. ^a). 

The finite-size scaling of D{N) can be better seen if 
this quantity is plotted as a function oi l/N rather than 
as a function of N. Such a plot is shown in Fig. ^h), 
including D{M, N) for fixed values of M. While the 
data points for M — 256 and M = 512 behave as if 
they would go to a finite value in the N ^ oo limit, the 
data available for M = 1024 (N < 256) do not show 
an upward curvature in 1/A^. The same is true for the 
extrapolated values of the dimer order parameter indi- 
cated by the solid line in Fig. EJb). Extrapolation us- 
ing Eq. ((ni) yields D{M ^ oo,iV -J> oo) = 0.003(2) 
and a very large correlation length, ^ — 5972. Since 
a flattening off is not apparent in the curves, we have 
repeated the fit without the exponential term, yielding 
D{M -> oo, oo) = 0.013(5). In order to obtain an 

upper bound for D, we have also carried out a fit to a 
second-order polynomial in l/N. The extrapolated value 
of D^hiM ^ oo, ^ oo) is 0.052(2) but this fit, indi- 
cated by the dashed line, clearly overestimates the bond 
order parameter significantly. 

By repeating the same analysis for several points in the 
vicinity of J± = 0.38 keeping Jx = 0.2 fixed we have also 
studied the finite-size scaling of D{N) in more detail and 



FIG. 3: (Color online) (a) Scaling of the dimer order param- 
eter as a function of the number of block states retained, M, 
for various chain lengths. The solid line is a fit using the 
function D{M,N) = D + a{N)M-l^'-^\ (b) The D{N) data 
of Fig. [2] plotted as a function of 1/A''. The solid line is a fit 
to the D{M oo,N) data using Eq. ITT)). The dotted line is 
a similar fit without the exponential factor, while the dashed 
line is a fit to a second-order polynomial in 1/7V. 



determined D{M oo, iV — j- oo). The results are shown 
in Fig. m For all data points extrapolation using Eq. pT|) 
yields D{M ^ oo, — > oo) ~ ©(lO^'') and a correlation 
length of the order 10^, except for = 0.38. It is ob- 
vious from the figure that the infiection point shows up 
only at larger system sizes as J±_ = 0.38 is approached, 
thus longer and longer systems must be investigated for 
a reliable extrapolation. It is worth mentioning, how- 
ever, that in all cases the extrapolated curves start with 
a horizontal slope. 

Since the decrease of the dimer order parameter is 
the slowest at J±^ = 0.38 we have reinvestigated this 
point using a better extrapolation scheme, i.e., the DBSS 
method. By setting a limit on the quantum information 
loss, the number of block states is selected dynamically 
during the RG steps of the DMRG. In our calculations, 
the maximum number of block states, Afmax, varies in 
the range 200 — 2500 for decreasing x- For x — lO^'', we 
have taken Mmin = 64, while we taken Mmin = 256 for 
all other cases. The largest dimension of the superblock 
Hamiltonian was over 14 million. The number of block 
states needed to achieve the predetermined x value is 
shown as a function of 1/A^ in Fig.[SJa). Obviously Mmax 
diverges for increasing N and decreasing x- For example, 
M ~ 2500 states must be kept for a ladder with TV = 256 
rungs in order to obtain an accuracy of x = 10^*". 

The dimer order parameter calculated for different ac- 
curacy levels is plotted as a function of x in FiglSfb) 
on a semi-logarithmic scale, while these data and the 
extrapolated values D(x ~^ 0, N) are plotted as a func- 
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oo) = 0.054(2). 
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FIG. 4: (Color online) Scaling of D{M — >■ oo, A*') as a function 
of for various J± values in the vicinity of J± — 0.38 
keeping Jx = 0.2 fixed. The solid line is a fit to the D{M -> 
00, A'') data using Eq. (fTTj). 
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FIG. 5: (Color online) (a) Maximum number of block states 
Mmax selected dynamically for given accuracy thresholds x 
plotted as a function of inverse system size 1/A'^. (b) Scaling 
of the dimer order parameter as a function of the quantum 
information loss x for various chain lengths. 



tion of 1/N in Fig. IH^a). Extrapolation to iV — cxd was 
again performed using Eq. yielding D{x — )■ 0, — J> 

oo) — 0.003(1), while the fit without the exponential 
term yielded D{x ^ 0,iV -J> oo) = 0.021(5). A portion 
of Fig. [nija) is displayed in an enlarged scale in Fig. ^h). 
It can be clearly seen that a small error, in this case of 
the order of 10"'', in D{x — >■ 0,N) can lead to a sig- 
nificant shift of the extrapolated value. We have again 
determined an upper bound through a fit to a second- 
order polynomial in 1/A^, obtaining D^hix ^ 0, N 
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FIG. 6: (Color online) (a) Similar to Fig. EJb) but with data 
calculated with the indicated fixed values of x and extrapo- 
lated to X = using an exponential form as shown in[5l^b). 
The sohd hue is a fit to the D{N, x — 0) data using Eq. fH} . 
The dotted line is a similar fit without the exponential factor, 
while the dashed line is a fit to a second-order polynomial in 
1/N. (b) Extrapolated data and fits for an enlarged region. 

To compare the results obtained from the M — >■ oo 
extrapolation with those of the x ^ extrapolation, we 
display both D{M -> oo,N) and D{x 0,N) on the 
same axes in Fig. [7]Ja). On the scale used, the results of 
the two procedures coincide. The difference is of order 
lO^'*, as can be seen on the logarithmic scale used in 
Fig.[7](b). The extrapolated data points show no upward 
curvature as a function of 1/N. Thus, we find that there 
is no evidence for a finite dimer order parameter to the 
resolution of the calculations performed here. 



V. CONCLUSION 

In summary, we have carried out density-matrix 
renormalization-group calculations on the s = 1/2 cross- 
coupled spin ladder, in which, in addition to antiferro- 
magnetic nearest-neighbor couplings along the legs and 
across the rungs of the ladder, an antiferromagnetic cou- 
pling between spins diagonally across the plaquettes is 
present. Our aim was to search for indications for the 
gapped columnar dimer phase found by Starykh and 
Balents,^^ in view of the fact that earlier numerical 
work22i^ only provided hints that this phase might ex- 
ist. Since it is known that the dimerized phase, if it 
exists, appears in a small parameter regime due to the 
marginally relevant character of the current-current cou- 
pling between the legs, we have performed the numerical 
calculations using the DMRG in a narrow region in the 
parameter space, namely for 0.36 < J± < 0.4, Jx = 0.2 
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FIG. 7: (Color online) (a) Similar to Fig.[3](b), but with the 
M — !> c« and x ~^ extrapolated data plotted on the same 
scale and (b) their difTerence plotted on a logarithmic scale as 
a function of 1/A'^. 



and have analyzed the accuracy of the calculation care- 
fully. Whereas the longest ladder studied in Ref. [lH had 
192 rungs and a maximum of 500 states were kept, we 
have carried out calculations for ladders with = 512 
rungs keeping M = 512 states and with N = 256 rungs 
keeping M = 1024 states. When the dimer order param- 
eter is plotted as a function of N (see Fig. [5]), the curves 
seem to tend to a finite value for large systems; the slope 
tends to zero. This was interpreted in Ref. 22, as an indi- 
cation of the existence of the dimerized phase. However, 
when more and more states are kept, this tendency be- 
comes less and less pronounced. When the dimer order 
parameter calculated for different M values is extrapo- 
lated to M cx), D curves downwards as a function of 
1/A^, as seen in Fig. ^h). In fact, finite-size scaling of 
the extrapolated D{M — > oo, N) data in the whole region 
0.36 <J±< 0.4, Jx = 0.2 using Eq. ^ yielded a dimer 
order parameter of order 10"'^ or smaller. 



In addition, at a particular point of the parameter 
space, J± = 0.38, Jx = 0.2, at which the phase is most 
likely to occur we have also performed calculations us- 
ing the DBSS method taking a number of different fixed 
values for the quantum information loss x for ladders of 
up to = 256 rungs keeping up to 2500 block states 
and performing an extrapolation to x = 0. One impor- 
tant message of this calculation is that if longer ladders 
are to be considered, the number of states kept in the 
DMRG procedure must be taken to be up to approxi- 
mately 10000 to obtain sufficient accuracy. This indi- 
cates the fundamental limitation of the DMRG method 
for the cross-coupled spin ladder. 

We have found that the result obtained for the dimer 
order parameter in the x = limit coincides with the one 
obtained from the M — > oo limit. The extrapolated data 
do not show upward curvature as a function of 1/N, in 
contrast to the results obtained on shorter chains with a 
smaller number of block states kept. The limiting value 
for iV — > oo depends strongly on the functional form used 
for the extrapolation. According to our best estimate, 
the maximum possible value of the dimer order parame- 
ter is of the order of the error of the calculation. While 
it cannot be completely ruled out that the curvature of 
D{N) as a function of 1 /N changes at much larger system 
sizes if there is a very small dimerization gap and a very 
large correlation length, there is no indication of such 
behavior in the present calculations. Therefore, we find 
that whether or not the columnar dimer phase is present 
anywhere in the phase diagram of the cross-coupled an- 
tiferromagnetic ladder remains an open question. 
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